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Qh! Abstract 



The .fsT — ► 7T7T decay amplitudes are studied within the framework of generalized factorization 
in which the effective Wilson coefficients are gauge- invariant, renormalization-scale and -scheme 



independent while factorization is applied to the tree-level hadronic matrix elements. Nonfactorized 
contributions to the hadronic matrix elements of (V — A){V — A) four-quark operators, which 
are needed to account for the suppression of the AI = 3/2 K — > tttt amplitude A2 and the 
^2 I enhancement of the AI = 1/2 ^o amplitude, are phenomenologically extracted from the measured 

K + — » 7T + 7r° decay and found to be large. The Aq/ 'A2 ratio is predicted to lie in the range 15-17 
for m s (lGeV) = (127 — 150) MeV. Vertex and penguin-type radiative corrections to the matrix 
elements of four-quark operators and nonfactorized effects due to soft-gluon exchange account for 
the bulk of the AI = 1/2 rule. Comparison of the present analysis with the chiral-loop approach 



is given. 



I. INTRODUCTION 

The effective Hamiltonian approach is the standard starting point for describing the 
nonleptonic weak decays of hadrons. In this approach, the decay amplitude has the form 
A ~ X) Cj(/x)(Oj(/i)), where the renormalization scale \x separates the short-distance contribu- 
tions contained in the Wilson coefficient functions Ci(fi) and the long-distance contributions 
contained in the hadronic matrix elements of 4-quark operators (Oi(fi)). Of course, the phys- 
ical amplitude should be independent of the choice of the renormalization scale and scheme. 
This means that the matrix elements have to be evaluated in the same renormalization 
scheme as that for Wilson coefficients and renormalized at the same scale fi. In principle, 
the scale fi can be arbitrary as long as it is large enough to allow for a perturbative calcula- 
tion of Wilson coefficients. In practice, it is more convenient to choose /i to be the scale of 
the hadron mass of the decaying particle so that the logarithmic term in the matrix element 
(O(fi)), which is of order ln(M 2 / fi 2 ) with M being the hadron mass, is as small as possible, 
leaving the large logarithms, which are summed to all orders via the renormalization group 
technique, to c(fi). 

Since the hadronic matrix elements are very difficult to calculate, it is not surprising that 
the issue of their fi dependence is generally not addressed in the literature. For meson decays, 
a popular approach is to evaluate the matrix elements under the factorization hypothesis so 
that (O(fi)) is factorized into the product of two matrix elements of single currents, governed 
by decay constants and form factors. However, the information of the scale dependence of 
(O(fi)) is lost in the factorization approximation because the vector or axial- vector current is 
partially conserved and hence scale independent. Consequently, the fi dependence of Wilson 
coefficients does not get compensation from the matrix elements. Although the correct fi 
dependence of (0(/i)) should be restored by the nonfactorized contributions to hadronic 
matrix elements, the difficulty is that nonfactorized effects are not amenable owing to their 
nonperturbative nature. Hence, the question is can we apply factorization to the matrix 
elements and in the meantime avoid the scale problem with (O(fi)) ? 

Fortunately, the fi dependence of hadronic matrix elements is calculable in perturbation 
theory. After extracting the fi dependence of (O(fi)) and combining it with the Wilson 
coefficients c(fi), we obtain renormalization-scale and -scheme independent effective Wilson 
coefficients. Then the factorization approximation can be safely applied afterwards to the 
matrix elements of the operator 0{ at the tree level. 

For the case of kaon decays such as K — > nir, it is obvious that fi cannot be chosen to 
be of order vhk] instead it has to be at the scale of 1 GeV or larger so that c(fi) can be 
reliably computed. Conventionally, the fi dependence of matrix elements involving kaons 
and pions are calculated by considering chiral loop corrections to (O). When chiral loops are 
regularized using the same dimensional regularization scheme as that for Wilson coefficients, 
the fi dependence of long-distance contributions will presumably match the scale dependence 
of Wilson coefficients so that the resulting physical amplitude is fi independent. While the 
scale dependence of K — > im matrix elements can be furnished by meson loops, it is clear that 
this approach based on chiral perturbation theory is not applicable to heavy meson decays. 
Therefore, it is desirable to consider the nonleptonic decays of kaons and heavy mesons within 



the same framework of generalized factorization in which the effective Wilson coefficients 
df are renormalization-scale and -scheme independent while factorization is applied to the 
tree-level hadronic matrix elements. The purpose of the present analysis is to see if our 
understanding of the AI = 1/2 rule can be improved in the effective Hamiltonian approach. 

The celebrated AI = 1/2 rule in kaon decays still remains an enigma after the first 
observation more than four decades ago. The tantalizing puzzle is the problem of how to 
enhance the Aq/A 2 ratio of the AI = 1/2 to AI = 3/2 K — > ttit amplitudes from the 
outrageously small value 0.9 [see Eq. ( |3.32|) below] to the observed value 22.2 ± 0.1 (for a 
review of the AI =1/2 rule, see JI[). In the approach of the effective weak Hamiltonian, 
the Aq/A 2 ratio is at most of order 7 even after the nonfactorized soft-gluon effects are 
included |L| . Moreover, the /x dependence of hadronic matrix elements is not addressed in the 
conventional calculation. In the past ten years or so, most efforts are devoted to computing 
the matrix elements to 0(p A ) in chiral expansion. This scenario has the advantages that 
chiral loops provide the necessary scale dependence for hadronic matrix elements and that 
meson loop contributions to the Aq amplitude are large enough to accommodate the data. 
However, it also becomes clear that chiral loops alone cannot explain the A 2 amplitude 
(see Sec. IV). Consequently, it is necessary to take into account nonfactorized effects on 
K + — > 7r + 7r° in order to have an additional suppression for the AI = 3/2 transition. 

Contrary to the chiral approach, the difficulty with the \x dependence of the physical K — ► 
mi amplitude is circumvented in the present analysis by working in the effective Hamiltonian 
approach in which the effective Wilson coefficients are gauge- invariant, renormalization- 
scale and -scheme independent. This approach is not only much simpler than chiral loop 
calculations but also applicable to heavy meson decays. By extracting nonfactorized effects 
from K + — > 7r + 7T°, we shall see that nonperturbative effects due to soft-gluon exchange 
and perturbative radiative corrections to four-quark operators account for the bulk of the 
observed AI = 1/2 amplitude. 

The present paper is organized as follows. In Sec. II we construct scheme and scale 
independent effective Wilson coefficients relevant to kaon decays. The K — > txtx matrix 
elements are evaluated in Sec. III. Comparison of the present analysis of the AI = 1/2 rule 
with the chiral loop approach is made in Sec. IV. Sec. V is for the conclusion. 

II. FRAMEWORK 

The effective Hamiltonian relevant to K — > txtx transition is 

C ( 10 \ 

H cS (AS= 1) = -^V ud V: s i 5>(aOO<M J +h.c, (2.1) 

where 

Ci (/•*) = Zi(fi) + Tyi(fj), (2.2) 

with r = -V td V*/(V ud V: s ), and 



0\ = (ud) v ^ A (su) v _ A , 2 = {u a bp) v _ A {q p u a ) v _ A) 

03(5) = (^)v-aS(w)v-a(v+x)> °4(6) = Wlv-AH^^V-il^+A)! ( 2 - 3 ) 
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°7(9) = «( 5f Ov-A H e g(w)v+A(v-A) ? °8(10) = ~ (Sad/?)v-A H e q(Wa) v+A ( v . A ), 

Z q Z q 

with O3-O6 being the QCD penguin operators, O7-O10 the electroweak penguin operators 
and (qiQ2) v±A = 9i7/x(l ± 75)92- The sum in Eq. fl2.3|) is over light flavors, g = u,d,s. It 
is obvious that only the Wilson coefficients z% are relevant to our purposes, as we are only 
interested in the CP-conserving part of K — > tttt transitions. 

In order to apply the factorization approximation to hadronic matrix elements, we need 
to compute the /x dependence of matrix elements arising from vertex and penguin-type 
radiative corrections to four-quark operators and combine it with the Wilson coefficients to 
form renormalization-scale and -scheme independent effective Wilson coefficient functions 
(for details, see P,U): 



4 s = z M + g (tT ^ l f + %)*M> 



^= Z M+^[ir^^+r T v j^), 



zf = zM + g (l ( r In & + r^ ZiM - ^(C t + C p ), 



1 \ a ^ 



zf = *M + r (lT In ^ + 4] Ziijx) +^(Q + C p ), 
4vr \ fi / 4 . 8tt 

*f = *(ri + £ (7? )T in f + f?) *o*) - ^(a + c p ; 

*f = *,(*«) + g (t^ In ^ + r^ zifr) + £& + C p ), 

zf^zM + ^^r^^ + r^z^), 

zf = zM + £ U )T In ^ + fi] ztM + ^C e , 

4vr \ fi J 9 . 8vr 

zf^zM + ^L^ln^ + rT) ztM, (2.4) 

where the superscript T denotes a transpose of the matrix, the anomalous dimension matrix 
7y as well as the constant matrix fy arise from the vertex corrections to the operators 
0\ — O10, C t , C p and C e from the QCD penguin-type diagrams of the operators Oi 2 , the 
QCD penguin-type diagrams of the operators 3 — 6 , and the electroweak penguin-type 
diagram of Q\ 2 , respectively: 



C t = G(m u )zi, 

C p = [G(m s )+G(m d ))z 3 + J2 G(m i )(z i + z 6 ), 

i=u,d,s 
O 

C e = -G(m u )(zi + 3z 2 ), 

2 
G{m q ) = -K-G(m q ,k,ii), (2.5) 

with k being a parameter characterizing the 75 scheme dependence in dimensional regular- 
ization, for example, 

in the naive dimensional regularization (NDR) and 't Hooft-Veltman (HV) schemes for 75. 
The function G(m, k, //) in Eq. (|2.5|) is given by 

G(m,k,fi) = -4 : dxx(l-x)\nl ^ '- , (2.7) 

where /c 2 is the momentum squared carried by the virtual gluon. 

The matrix f in (|2.4| ) gives momentum-independent constant terms which depend on the 
treatment of 75 in dimensional regularization. An early evaluation of f is performed in the 
off-shell quark scheme ||. However, it was pointed out by Buras and Silvestrini that 
zf thus constructed suffer from gauge and infrared ambiguities since an off-shell external 
quark momentum, which is usually chosen to regulate the infrared divergence occurred in the 
radiative corrections to the local 4-quark operators, will introduce a gauge dependence. It was 
shown recently in that the above-mentioned problems on gauge dependence and infrared 
singularity connected with the effective Wilson coefficients can be resolved by perturbative 
QCD (PQCD) factorization theorem. In this formalism, partons, i.e., external quarks, are 
assumed to be on shell, and both ultraviolet and infrared divergences in radiative corrections 
are isolated using the dimensional regularization. Because external quarks are on shell, gauge 
invariance of the decay amplitude is maintained under radiative corrections to all orders. This 
statement is confirmed by an explicit one-loop calculation in [[|. The obtained ultraviolet 
poles are subtracted in a renormalization scheme, while the infrared poles are absorbed 
into universal nonperturbative bound-state wave functions. Explicitly, the effective Wilson 
coefficient has the generic expression 

c cS = c(fi) gi (fi)g 2 (fi f ) , (2.8) 

where #i(/i) is an evolution factor from the scale /i to rriQ, whose anomalous dimension 
is the same as that of c(fi), and ^(/^z) describes the evolution from rriq to \if (/i/ being 
a factorization scale arising from the dimensional regularization of infrared divergences), 
whose anomalous dimension differs from that of c(/i) because of the inclusion of the dynamics 
associated with spectator quarks. For kaon decays under consideration, there is no any heavy 
quark mass scale between m c and txik- Hence, the logarithmic term emerged in the vertex 



corrections to 4-quark operators is of the form In////// as shown in Eq. ( |2.4|) . We will set 
fif = 1 GeV in order to have a reliable estimate of perturbative effects on effective Wilson 
coefficients. 

The scale dependence of vertex and penguin-type corrections shown in Eq. ( |2.4| ) is 
governed by the terms 7^ In// and G(/i), while the 75-scheme dependence is determined by 
the matrix fy as well as G(k). Formally, one can show that the // and 75-scheme dependence 
of the next-to-leading order (NLO) Wilson coefficient, say zi(fi), is compensated by 7^ In// 
and fy, respectively, to the order of a s [§,[§. This means that the NLO Wilson coefficients 
Zj(/i) appearing in Eq. ( |2.4j ) together with a s /4ii or a/47r should be replaced by the lowest- 
order values zf (fi). The numerical values of z? are displayed in Table I. We see that 
except for z$ 7 , effective Wilson coefficients shown in the last two columns of Table I are 
indeed renormalization scheme independent, as it should be. 

TABLE I. AS = 1 Wilson coefficients at u = 1 GeV for m t = 170 GeV and A^ = 325 

^ l MS 

MeV, taken from Table XVIII of [6] . Also shown are the effective Wilson coefficients obtained from 
zf DR {fx) and zf Y {{i) via Eq. (2.4) with fi = 1 GeV, /// = 1 GeV and k 2 = m 2 K /2. 





LO 


NDR 


HV 


zf(NDR) 


zf(HV) 


Z\ 


1.433 


1.278 


1.371 


1.718 


1.713 


Z-2 


-0.748 


-0.509 


-0.640 


-1.113 


-1.110 


ZZ 


0.004 


0.013 


0.007 


0.034 


0.033 


ZA 


-0.012 


-0.035 


-0.017 


-0.088 


-0.087 


Zf> 


0.004 


0.008 


0.004 


0.026 


0.026 


z% 


-0.013 


-0.035 


-0.014 


-0.093 


-0.089 


z^/a 


0.008 


0.011 


-0.002 


0.063 


0.069 


z 8 /a 


0.001 


0.014 


0.010 


0.016 


0.013 


zg/a 


0.008 


0.018 


0.005 


0.072 


0.078 


zio/a 


-0.001 


-0.008 


-0.010 


-0.011 


-0.012 



In the late 70's and early 80's, it had been suggested that penguin operators may account 
for the AI = 1/2 rule observed in kaon decays. With the advent of the effective Hamiltonian 
approach, it is realized that the AI = 1/2 selection rule cannot be dominated by the penguin 
mechanism. One popular argument is that at the scale, say // > m c , the penguin Wilson 
coefficients are negligible due to the (incomplete) GIM mechanism. Since the penguin Wilson 
coefficients become important at // = 1 GeV, for instance, one may wonder if the physical 
penguin contributions to K — > nn is independent of the choice of //. The point is that 
although z 3 , ■ ■ ■ , zio vanish at, say /i = 2 GeV, the effect of the penguin diagrams with the 
internal u quark induced by the current-current operator 0\ has to be taken into account 
when evaluating matrix elements. Consequently, the total penguin contribution is scale 
independent, and this is the merit of effective Wilson coefficients in which the perturbative 
effect of the penguin diagram with the internal u quark is already included. 
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We add a remark before ending this section. A dynamical phase can arise from the 
time-like penguin diagram involving internal u loop quarks. However, since this penguin- 
induced phase is incorporated into the isospin zero final-state interaction phase shift So to 
be introduced below in Eq. ( |3.11| ), its contribution should not be double-counted in the 
effective Wilson coefficients. 



III. CALCULATIONS 

In this section we first study the K —* tttt matrix elements based on the vacuum insertion 
approximation, and then turn to the AI = 1/2 and AI = 3/2 amplitudes. 

A. Matrix Elements 

It is convenient to make isospin decomposition of the matrix elements (Oj)o,2 = (tttt, I = 
0, 2\Oi\K°) which are related to K — tttt transitions via 

(0<) O = 4= (2(TT + TT-\O t \K°) + (TT°TT°\O t \K ) 

(o z ) 2 = -L ((tt+tt-io.ik ) - (ttV|o<|jr°>) = y|(7T + 7r°|o^ + ). (3.1) 

Conventionally, the matrix elements (Oi)o,2 are evaluated using the vacuum insertion ap- 
proximation (i.e. factorization hypothesis). Under this assumption, we have, for example, 

(7T + 7r-|Oi|^°) = (TT + \(ud)\0}{TT-\(su)\K ) + -^ (tt+tt" | (uu) |0) <0| (sd) \K°) 

= fAm 2 K -ml)F ^(ml), * (3.2) 

where the form factor Fq is defined in and the ^-exchange contribution vanishes due to 
vector current conservation. The q 2 dependence of the form factor Fq is usually assumed to 
be dominated by near poles in a monopole manner: 

Fttf) = ^^ - i?*(0) (l + £) , (3-3) 

where m* is the pole mass of the + scalar meson with the quantum number of sq (q = u,d). 
In chiral perturbation theory (ChPT), we have F^ n (0) = 1 due to vector current conservation 
and (see [I] for details) 



m 



^ = ^=f^-l)^^»^ (M 



fl \U J (m 2 K - ml) AI 



where L5 is one of the coupling constants in the 0(p 4 ) chiral Lagrangian for strong inter- 
actions, and A x p» 2nfn (our /„. = 132 MeV) is the chiral-symmetry breaking scale ||. 
Therefore, 
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(7T + 7r-|0 1 |^°)= U(m 2 K -ml)[l 



ml 



(3.5) 



In ChPT, the term proportional to m 2 /A 2 is counted as a contribution of 0(p 4 ). 

Contrary to charmless B decays, the penguin operators 56 , do not contribute directly 
to K° — > 7r7r because of the wave function 7r° = -j^[uu — dd) and the SU(3)-singlet nature 
of the V + A current. Nevertheless, the Fierz transformation of O^q via (V — A)(V + A) — > 
—2(S + P)(S — P) does make contributions. For example, 



(vr + 7i- \0 6 \K°) = -(7r + |M7 5 d|0)(7r-|sw|ir°) - -(7r + 7r-|dd|0)(0|s7 5 ci|X ) 
o o 

+ ^({dd) + (ss))(n+n-\s l5 d\K°). 



(3.6) 



The matrix elements of scalar and pseudoscalar densities can be evaluated using equations 
of motion or the chiral representation of quark densities. The former method gives 



(n-(q)\su\K°(k)) = v 



1 + 



(n+(q + )n-(q_)\dd\0)= v 



1 + 



(k - q)' 

A 2 
x 

(g+ + q-f 

A x 



(tt + (q)\wy 5 d\0) = if w v, 
(0\s l5 d\K°(k))=tf K v, 



(3.7) 



where uses of Eqs. ( |3.3| ) and ( |3.4| ) have been made, and 



2 



m 2 K0 



m u + m d m d + m s 



m K- m * 

m s — m„ 



(3.8) 



characterizes the quark-order parameter (qq) which breaks chiral symmetry spontaneously. 
The second term on the r.h.s. of Eq. ( |3.6| ) is the so-called spacelike penguin contribution. 
Unlike the case of hadronic charmless B decays, the spacelike penguin diagram in K — > tttt 
is calculable. The last term in Eq. (|3.6| ), which is a tadpole contribution arising from the 
vacuum expectational values of quark bilinears, does not contribute to the physical K — > tttt 
amplitude M. Hence, we obtain 



(n + (q + )n-(q_)\0 6 \K°(k)) = -Qf^v 2 *" ' r ~ 



4 
',— 
3* 



A x 



o 



M 



(3.9) 



where we have applied Eq. 

The matrix elements obtained under the vacuum insertion approximation are summarized 
below: 



(O 2 ) 
(O 3 )o 



1 



-X 



N r 



N r . 



-X, 



(o 1 ) 2 = 


3 


(- 


1 


(o 2 ) 2 = 


-^x 

3 


(> + 


1 


(O 4 )o = 


--x, 







io^fj^^-x, { o 7h = fjy-±=x, 

(o 8 )o = VQf K v 2 + ^-x, <o 8 > 2 = Vsfy - j^j=X, 

(09)0 = -\x (i-±), (o 9 ) 2 = -±=x (1 J 



(o u >„ = ix(i- ] L) > (o 10 ) 2 ^-Lx(i + ^), (3.10) 



where X = v/3/2 firiw^ — m^) and 1/A£ corrections to (V — A)(V ± A) matrix elements as 
well as 1/A^ corrections to (S* + P)(5 — P) matrix elements have been neglected. 

In terms of the isospin matrix elements, the corresponding isospin decay amplitudes are 
given by 

A ie * 5 > = ^=V ud V: s J2c^)(0^))r, 

= ^v ud v: s J24 a (o l ) I , (3.11) 

and hence 

ReA , 2 = ^^5>f<O,,> , 2 , (3.12) 

V2 cosd , 2 i 

where 5o and 5 2 are S-wave im scattering isospin phase shifts. In the present paper, we will 
use the analysis of |§ for phase shifts: 

5 = (34.2 ±2.2)°, S 2 = -(6.9 ±0.2)°. (3.13) 

Experimentally, the isospin K — > tttt amplitudes are given by 

Re A = 3.323 x 1(T 7 GeV, Re A 2 = 1.497 x 10" 8 GeV. (3.14) 

B. The AI = 3/2 amplitude 



It is straightforward to show from Eqs. ( pM0| ) and ( |3.12| ) that the isospin 2 amplitude of 
K — > 7T7T has the form 



.4 



(o) 



G F V ud V*[\ ,3, xl y/2 



y/2 cos5 2 
where 



ai + a 2 + -(—07 + a 9 + a Wj 



3 -X + V3Uv 2 a 8 , (3.15) 



a 2i - ^ 2 j + jy^^j-1, «2J-1 - ^ 2 j_l + ]y~ Z 2J (3.16) 



for i = 1, • • • , 5, and the superscript (0) indicates that this amplitude is induced by pure 
AI = 3/2 weak interactions. Since the electroweak penguin coefficients are very small 
compared to zf and zf , it is clear that the AI = 3/2 decay amplitude is entirely governed 
by current-current 4-quark operators. It is also known that K + — > 7r + 7r° (or AI = 3/2 
K° — > 7T7T decays) can be generated from the AI =1/2 decays K + — > 7i + r] (r/) followed by 
the isospin breaking mixing tt° — r] — rj' [|H],|Il]]. As a result, the total AI = 3/2 amplitude 
reads 



A, 



4 0) 



(3.17) 



where the expression of Qm = A 2 B /A 2 can be found in Appendix A . Employing the quark 
mass ratios m^jm u = 0.553 ± 0.043 and ms/rrid = 18.9 ± 0.8 obtained in a recent detailed 
analysis based on ChPT |I2 |, we find from Eq. ( |A7| ) that 



n 



IB 



0.25 ±0.02. 



Eqs. (|3~T5|)-(|3~T8|) lead to 

A 2 = 4.133 (zf + z. 



eflf\ 

2 ) 



1 + — ) x l(r s (M\' 

■L V r. 



(3.18) 



(3.19) 



Using the effective Wilson coefficients zf E given in Table I, it is easily seen that the predicted 
A 2 is too large by a factor of 2.2 compared to experiment (|3.14|) . This means that nonfac- 
torized contributions that have been neglected thus far should be taken into account. For 
K — > 7T7T decays, nonfactorizable effects in hadronic matrix elements can be absorbed into 



the parameters a| [IT3HTH]: 



hi 



Z 2i + ( jr + X2i ) Z2i-1 ■ 



7 Cff 

hi-i 



Z 2i-l + ( J7 + X2i~l ) Ai 1 



(3.20) 



where the nonfactorized terms Xi,2 relevant to K + —* -k + tt° decay are given by [|T~5| [3~7 1 

_(_ft:+7r+,7r ) a 2 (K+TT+,n°) 



{K + TT°,7T+) Ol (K+TT a ,7T+) 

z 2 



X2 



z l 



(3.21) 



with ai : 2 = zl * + zf[/N c , and 

{K+iP,*+) __ (K + K°\(ud) v _ A (su) v _ A \K + } nf _ (7T + 7r°\(ud) v _ A (su) v _ A \K + ) 



(Tr+ir°\(ud) v _ A (su) v _ A \K+) f (n+\(ud) v _ A \0}(^\(su) v _ A \K+) 
+K+*°,*+) __ I (7r+7r°\(uX a d) v _ A (sX a u) v _ A \K+) 



-1, 



(3.22) 



2 (ir+\(ud) v _ A \0)(ir°\(su) v _ A \K+)> 
being nonfactorizable terms originated from color-singlet and color-octet currents, respec- 



tively, and (qiX a q 2 ) v _ A = giA a 7 M (l — 75)^2- The subscripts 'f and 'nf in Eq. ( |3.22| ) stand for 
factorizable and nonfactorizable contributions, respectively, and the superscript (K + tt°, tc + ) 
in Eq. ( |3.21 ) means that the ir + is factored out in the factorizable amplitude of K + — > 7r + 7r° 
and likewise for the superscript (K + ir + , n°). In the large- N c limit, E\ = 0(1/ N^) and 
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^8 = 0(1/N C ) |17|. Therefore, the nonfactorizable term x in the iV c — »• 00 limit is dominated 
by color octet-octet operators. 

Assuming xi — X2 in Eq. (|3.20| ) for af and af and fitting Eq. ( [3.19|) , in which l/N c is 
replaced by 1/N C + x, to the experimental value ( |3.14| ), we obtainF] 



X (K -> tttt) = -0.73 , (3.23) 

and hence 

af = 2.16, af = -1.80. (3.24) 

For comparison, the nonfactorized effects in hadronic two-body decays of charmed and bot- 
tom mesons are given by |I8" | 



X 2 (D -»• Kir) ~ -0.33 , x 2 (-B -► -D*") ~ (0.12 - 0.21) . (3.25) 

The fact that 

\ X (K -► tttt)! » \ X2 (D -► Ftt)| » 1x2(5 - Dtt)| (3.26) 

is consistent with the intuitive picture that soft gluon effects become stronger when final- 
state particles move slower, allowing more time for significant final-state interactions after 



hadronization [13]. 



Note that in B or D decays, the parameters af 2 and hence xi,2 in principle can 
be determined separately from experiments under some plausible assumptions. For ex- 
ample, Xi{D — > Kit) and Xz{D — > Kit) can be extracted from the isospin analysis of 
D° — ► K~tt + , K tt° and D + —>■ K tt + data provided that the H^-exchange is negligible. [] 
By contrast, Xi{K ~^ 7r7r ) an d XiiK — > tttt) cannot be determined from the data without 
invoking a further assumption because neutral K° — > tttt decays receive additional pen- 
guin contributions. That is why we make the universality assumption xi — X2 to extract 
01,2 (-K" ~^ tttt) from the measurement of K + — > tt + tt° . 

In the literature, the effective parameters af are sometimes expressed in terms of the 
the scheme- and scale-dependent Wilson coefficients Zj(/i), for example, 

af = zM + (— + xM) Mfi), of = Mfi) + (jf + X2O*)) zM, (3.27) 

where we have put a tilde on xi,2 to distinguish them from xi,2 defined in Eq. ( |3.20| ). Then 
it is clear that xi,2 must be 7 5 -scheme and scale dependent in order to ensure the scheme 
and scale independence of af 2- Comparing Eqs. (|3.27|) and ( |3.20| ), we see that xi,2 receive 
contributions from vertex radiative corrections. It should be stressed that the assumption 



*In the so-called large- N c approach, one has x = ~ 1/3 to the leading 1/N C expansion, 
t Since in general \z\/z2\ 3> 1, the determination of X2 is easier and more reliable than xi- 
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Xi — X.2 cannot lead to 75-scheme independent af|. To see this, we assume ( |3.24|) to be the 



true values for a\ s 2 an d apply the Wilson coefficients evaluated at fi — 1 GeV in NDR and 
HV schemes given in Table I. We find 

xrV) = -2-07, ^ DR (/i) = -l-34, 

X? V G") = -l-57, ^ V (/i) = -l-18, (3.28) 

at /i = 1 GeV by fitting ( p. 27] ) to ( |3.24| ). This implies that phenomenologically it is not 
possible to determine a\ 2 from the data of K — > titt if we start with the scheme- and 
scale-dependent Wilson coefficients Zj(/i) without taking into account vertex corrections to 

Xl,2- 

C. The AI = 1/2 amplitude 



From Eqs. ( |3.10| ) and ( |3 . 1 2| ) we obtain the AI =1/2 amplitude: 



An = ^ 



g f v ud v: 



211. . 

-ai - -a 2 + a 4 + -{a 7 -a 9 + a 10 ) 



X 



\/2 cos5 
- 2V6v 2 (f K - Qa e + V6v 2 f K aX (3-29) 

where we have neglected the contribution arising from it — 77 — 7/ mixing. For simplicity, we 
have also dropped the superscript 'eff' of the parameters a«. To incorporate nonfactorized 
effects, we shall make the universality assumption: 

Xll = Xi = X2 = X3 = X4 = X9 = X10, 

Xlr = X5 = Xg = X7 = Xs- (3.30) 

The nonfactorized effects in the matrix elements of (V — A)(V + A) operators are a priori 
different from that of (V — A)(V — A) operators. Indeed, we have learned from hadronic 
charmless B decays that xlr 7^ Xll @- However, in the absence of information for the 
nonfactorized contributions to K — » nn penguin operators, we shall assume xlr ~ Xll = 
—0.73 for simplicity. Moreover, we found in actual calculations, A is insensitive to the 
value of Xlr- From Eq. ( |3.29| ) and Table I, it is easily seen that the nonfactorized term 
Xll = —0.73, which is needed to suppress A 2 to the observed value, will enhance the tree 
contribution to A by a factor of 1.9; that is, the tree contribution to A /A 2 ratio is increased 
by a factor of 3 ! 

Treating the strange quark mass m s and hence the parameter v as a free parameter, we 
plot in Fig. 1 the ratio A /A 2 as a function of m s at the renormalization scale \x = 1 GeV. 
Specifically, we obtain 

Ao f 17.1 at m s (lGeV) = 127 MeV, 

A~2~\ 15.3 at m s (1 GeV) = 150 MeV. ^ ' ' 

It is clear that m s is favored to be smaller. Presently there is no consensus regarding the 
values of light quark masses. It is interesting to note that several recent lattice calculations 
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0.08 0.1 0.12 0.14 0.16 0.18 0.2 

m s ( GeV ) 

FIG. 1. The ratio of Aq/A 2 versus m s (in units of GeV) at the renormalization scale [i = 1 
GeV. The solid thick line is the experimental value for Aq/A 2 . 

give a lighter strange quark mass: Results using the Sheikholeslami-Wolhlert fermion yield 
m s = (95 ± 16) MeV [^9|1 , a computation based on domain wall fermions obtains m s = 
(95 ± 26) MeV pOfl , a quenched QCD calculation together with the quark mass ratios from 
ChPT gives m s = (97 ± 7) MeV J2TJ, and a new unquenched lattice result indicates a still 
lower number m s = (84 ± 7) MeV ||22|| , all in the MS scheme at a scale of 2 GeV. The strange 
quark mass 95 MeV at // = 1 GeV corresponds to m s = 127 MeV at // = 1 GeV. 

It is instructive to see how the prediction of the A/ = 1/2 rule progresses at various 
steps. In the absence of QCD corrections, we have a 2 = |ai and a 3 = a 4 • • • = aio = under 
the vacuum insertion approximation. It follows from Eqs. ( pM5p and ( |3.29| ) that |T]| 



A 2 



4V2 



0.9 (in absence of QCD corrections). 



(3.32) 



With the inclusion of lowest-order short- distance QCD corrections to the Wilson coefficients 
Z\ and z 2 evaluated at // = 1 GeV, A /A 2 is enhanced from the value of 0.9 to 2.0, and 
it becomes 2.4 for m s (lGeV) = 127 MeV when QCD and electroweak penguin effects are 
included. This ratio is suppressed to 1.8 with the inclusion of the isospin-breaking effect, 
but it is increased again to the value of 2.1 in the presence of final-state interactions with 
5 Q = 34.2° and 5 2 = —6.9°. Replacing c^°(/i) by the effective Wilson coefficients df- , or 
equivalently replacing the LO Wilson coefficients by the NLO ones and including vertex and 
penguin-type corrections to four-quark operators, we find A /A 2 = 4.8 . Finally, the inclusion 
of nonfactorized effects on hadronic matrix elements will enhance A /A 2 to the value of 17.1 . 
In short, the enhancement of the ratio A Q /A 2 is due to the cumulative effects of the short- 
distance Wilson coefficients, penguin operators, final-state interactions, nonfactorized effects 
due to soft-gluon exchange, and radiative corrections to the matrix elements of four-quark 
operators. Among them, the last two effects, which are usually not addressed in previous 
studies (in particular, the last one), play an essential role for explaining the bulk of the 
A/ = 1/2 rule. In present calculations, penguin operators account for 35% of the A/ = 1/2 
rule for m s {\ GeV) = 127 MeV. 

Note that thus far we have neglected the PF-exchange effect, vanishing in he vacuum inser- 
tion approximation. Since the VF-exchange amplitude in charmed meson decay is comparable 
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to the internal PU-emission one ||23|| , it is conceivable that in kaon physics the long-distance 
contribution to W -exchange is as important as the external P^-emission amplitude. There- 
fore, the jy-exchange mechanism could provide an additional important enhancement of the 
A /A 2 ratio. 

IV. COMPARISON WITH CHIRAL APPROACH 

Since the scale dependence of hadronic matrix elements is lost in the factorization ap- 
proach, it has been advocated that a physical cutoff A c , which is introduced to regularize 
the quadratic (and logarithmic) divergence of the long-distance chiral loop corrections to 
K —* 7T7T amplitudes, can be identified with the renormalization scale fi of the Wilson coeffi- 
cients |p4fl . A most recent calculation along this line which includes 0{p i ) tree contributions 



indicates that while the isospin amplitude A is largely enhanced, the amplitude A 2 is 
highly unstable relative to the cutoff scale A c and it even changes sign at A c > 650 MeV 
|25|,26]. The large uncertainty for A 2 arises from the fact that the two numerically leading 



terms, the tree level and the one-loop quadratically divergent term, have approximately the 
same size but opposite sign. 

Since the scale dependence of Wilson coefficients is of the logarithmic type, it seems 
quite unnatural to match the quadratic cutoff with the \x dependence of c(fi). Therefore, it 
is necessary to use the dimensional regularization to regularize the chiral loop divergences 
and apply the same renormalization scheme, say the MS scheme, in order to consistently 
match the scale dependence of Wilson coefficients evaluated using the same regularization 
scheme. In this case, the inclusion of chiral loops will make a large enhancement for Aq and 
a small enhancement for A 2 p7j. However, as stressed in passing, the naive prediction of A 2 



in the absence of nonfactorizable effects is too large (by a factor 1.6 in our case) compared 
to experiment. Therefore, chiral-loop corrections to A 2 will make the discrepancy between 
theory and experiment even worse. Evidently, this indicates that not all the long-distance 
nonfactorized contributions to hadronic matrix elements are fully accounted for by chiral 



loops. Several authors H28|j26| have considered different models, for instance the chiral quark 
model or the Nambu-Jona-Lasinio model, to incorporate nonfactorized contributions arising 
from soft gluonic corrections. For example, the nonfactorized gluonic corrections computed 
in the chiral quark model amount to replacing 1/N C in the matrix elements (01-3)0,2 [see 
Eq. (|3T0|)] by H 



1 1 , , 1 ( N c A^(a s GG/n) \ 

parametrized in terms of the gluon condensate (a s GG/n). It is clear that the correction 
— ^<gg>/3 plays the same role as the nonfactorized terms Xi defined in Eq. Q3.20Q . This 
nonfactorized effect is important since it can suppress the A 2 amplitude. Using (a s GG/ir) = 
(334 ± 4 MeV) 4 f2"9fl , one obtains 5(gg) — 1.51. It is clear that the soft gluon correction, 



corresponding to Xi = X2 = —0.50, is large enough to revert the sign of the 1/N C term and 
thus suppress the A 2 amplitude. 
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In the present analysis, the nonfactorized contribution to the matrix elements of (V — 
A)(V — A) operators characterized by the nonfactorized terms xi,2 = —0.73 is comparable 
to that obtained in the chiral quark model. Therefore, in a rough sense, vertex and penguin- 
type radiative corrections to K — > im matrix elements in the effective Hamiltonian approach 
corresponds to chiral-loop contributions in the aforementioned chiral approach. However, 
the penguin contribution to the Aq amplitude in the latter approach is usually smaller than 
that in the former. For example, the 0§ operator contribution to Aq is about 20% in [2S|, 
and it is even smaller in other chiral-loop calculations. 

V. CONCLUSION 

We have studied K — > tttt decays within the framework of generalized factorization in 
which the effective Wilson coefficients are renormalization-scale and -scheme independent 
while factorization is applied to the tree-level hadronic matrix elements. Nonfactorizable 
contributions to the hadronic matrix elements of (V — A)(V — A) four-quark operators 
are extracted from the measured K + — > 7r + 7r° decay to be xx,2 = -0.73 which explains 
the suppression of the A/ = 3/2 K — > im amplitude A 2 and the enhancement of the 
A/ = 1/2 Aq amplitude. The A/ = 1/2 rule arises from the cumulative effects of the 
short- distance Wilson coefficients, penguin operators, final-state interactions, nonfactorized 
effects due to soft-gluon exchange, and radiative corrections to the matrix elements of four- 
quark operators. In particular, the last two effects are the main ingredients for the large 
enhancement of A with respect to A 2 . The A /A 2 ratio is predicted to lie in the range 15-17 
for m s (l GeV) = (127 — 150) MeV. Comparison of the present analysis of the A/ = 1/2 rule 
with the chiral-loop approach is given. 
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APPENDIX 



A. ISOSPIN BREAKING EFFECTS ON A/ = 3/2 K -> tttt AMPLITUDE 

In this Appendix we give an updated estimate of isospin breaking contribution to K^ 
7T + 7r° due to the n — rj — rf mixing. Writing 



a 2 = 4 0) + 4 B , 



(Al) 



we have jITJ 






1 A m d -m u 

A 2 ~ 3^ A 2 



m. 



(cos 9 — v2 sin 9) (cos 9 — v2 sin i 



+ (sin0 + V2cos0)(sin0H-V2 — *--cos6)—l 



P 



1 2 ~\ 

mi — m 



1 + 5 ~ ' ml, — ml 



(A2) 



where 9 is the r\ — rf mixing angle and the parameters p and 8, defined by 

( V s\H w \K°) = y|(l + 8)(-k°\H w \K°} 7 



( Vo \H w \K°) = -2J-p(n \H w \K ), 



(A3) 



measure the breakdown of nonet symmetry in K° — 7] transition and of SU(3)-flavor sym- 



metry in K° — 7] S} respectively. We can use the radiative decays K L — > 77 and it c 
constrain p and 8: 



A{K L -> 77) 
A(ir° — > 77) 



_4. 



m 



K 



m 



K 



_98_ 



c, 



77 to 



(A4) 



where 

C = i 



2 2 

m A K — m^ 
mi- — ml 



m\ - m%, 



l ~(l + 5)cos9 + 2J~psmH 

-(1 + <5) sin fl-2W-p cost 



/1 u 



4f - e 



lh sin9+2 nf cos9 



(A5) 



and g§ = 0.26 x 10 5 m 2 K is the coupling constant in the AS = 1 effective chiral Lagrangian. 
From the data of Kl — > 77 and tt° — > 77, we find |C| — 0.87. Using 

£ = -15.4°, f 8 /f n = 1.26, f /U = 1.17 (A6) 

determined phenomenologically |30| and 8 = 0.17 | 3~]| , we obtain p = 0.96. Numerically, we 
find that Q\b is almost insensitive to the values of 8, p and 9 as long as they are constrained 
by (. Hence, to a very good approximation, we obtain 

m d -m u 



n m = (10.45 ±0.05) 



m. 



(A7) 



where the experimental value of Aq/A 2 = 22.2 has been used. 



B. ANOMALOUS DIMENSIONAL AND CONSTANT MATRICES 

For reader's convenience, we list here the anomalous dimensional matrix "f v and the 
constant matrix fy appearing in Eq. 
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in the HV scheme. Note that the 66 and 88 entries of fy given in || are erroneous and have 
been corrected here. 
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